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REPRESENTATION OE SOLUTIONS AND LARGE-TIME 
BEHAVIOR FOR FULLY NONLOCAL DIFFUSION EQUATIONS 

JUKKA KEMPPAINEN, JUHANA SILJANDER AND RICO ZACHER 


Abstract. We study the Cauchy problem for a nonlocal heat equation, which 
is of fractional order both in space and time. We prove four main theorems: 

(i) a representation formula for classical solutions, 

(ii) a quantitative decay rate at which the solution tends to the fundamental 
solution, 

(iii) optimal L^-decay of mild solutions in all dimensions, 

(iv) L^-decay of weak solutions via energy methods. 

The first result relies on a delicate analysis of the definition of classical 
solutions. After proving the representation formula we carefully analyze the 
integral representation to obtain the quantitative decay rates of (ii). 

Next we use Fourier analysis techniques to obtain the optimal decay rate for 
mild solutions. Here we encounter the critical dimension phenomenon where 
the decay rate attains the decay rate of that in a bounded domain for large 
enough dimensions. Consequently, the decay rate does not anymore improve 
when the dimension increases. The theory is markedly different from that of 
the standard caloric functions and this substantially complicates the analysis. 

Finally, we use energy estimates and a comparison principle to prove a 
quantitative decay rate for weak solutions defined via a variational formulation. 
Our main idea is to show that the L^-norm is actually a subsolution to a purely 
time-fractional problem which allows us to use the known theory to obtain the 
result. 


1. Introduction 

We study the Cauchy problem for the diffusion equation 

d^{u{t,x)—u^)-\-Cu{t,x) = f{t,x) in IR+x 0 < a < 1, (1.1) 

where tio(x) = u(0,a:) is the initial condition, 9" denotes the Riemann-Liouville 
fractional derivative if a C (0,1) and £ is a nonlocal elliptic operator of order 
P C (0,2]. A standard example is the fractional Laplacian £ = (—A)' 2 . The 
equation is nonlocal both in space and time and we call such a parabolic equation 
a fully nonlocal diffusion equation. 

Our emphasis is on the decay properties, and for the space-fractional heat diffusion 
such questions have been studied, for instance, by Chasseigne, Chaves and Rossi 
in |13] as well as by Ignat and Rossi in |26) . For a more comprehensive account 
of the asymptotic theory in case a = 1, we refer to m- The decay of solutions 
and behavior of the Barenblatt solution for the space-fractional porous medium 
equation has, in turn, been studied by Vazquez in m- In the present paper, we 
extend these developments - concerning the fundamental solutions, representation 
formulas and decay properties - to the above fully nonlocal equation. For the case 
/3 = 2, see earlier works by Vergara and Zacher in [43] and by Vergara and the 
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present authors in [^. For the regularity theory of nonlocal equations in case 
a = 1 or /? = 2, we refer to [miziisiiiisiEninisiiis] and the references therein. 

Nonlocal PDF models arise directly, and naturally, from applications. Time frac¬ 
tional diffusion equations are closely related to a class of Montroll-Weiss continuous 
time random walk (CTRW) models and have become one of the standard physics 
approaches to model anomalous diffusion processes [HI HH HU 133]. For a detailed 
derivation of these equations from physics principles and for further applications 
of such models we refer to the expository review article of Metzler and Klafter 
in [34]. The fractional Laplacian arises in the modelling of jump processes and also 
in quantitative finance as a model for pricing American options [TBl HU] . The fully 
nonlocal diffusion equation, in particular, has been used in diffusion models, for 
instance, in m and m- 

Despite their importance for applications, the mathematical study of fully nonlocal 
diffusion problems of type (ED is relatively young. In a very recent paper Allen, 
Caffarelli and Vasseur [T] have studied the regularity of weak solutions to such 
problems. Even more recently, simultaneously to our work, Kim and Lim m 
have considered the behavior of fundamental solutions, whereas Cheng, Li and 
Yamamoto [14] have studied other aspects of the asymptotic theory. Apart from 
these papers, the study of the parabolic problem has mostly concentrated on the 
aforementioned cases a = 1 or /3 = 2. 

We point out that the nonlocal in time term in ED, with i9“ being the Riemann- 
Liouville fractional derivation operator, coincides (for sufficiently smooth u) with 
the Caputo fractional derivative of u, see ED below. The formulation with 
Riemann-Liouville fractional derivative has the advantage that a priori less reg¬ 
ularity is required on u to define the nonlocal operator. In particular, our formu¬ 
lation is exactly the one which naturally arises from physics applications, see for 
instance [31 equation (40)]. 

Our first main result considers a representation formula for classical solutions of the 
Cauchy problem for equation ED with C = (—A) 2 . In the process, we calculate 
the exact behavior of the fundamental solutions. 

Next, we show that the mild solutions, which are defined through the representation 
formula whenever its integrals are finite, tend to the fundamental solutions Z and Y 
- corresponding to the initial and forcing data, respectively - in with quantitative 
decay rates. Such results are nontrivial already for standard caloric functions, 
especially in the case of a non-vanishing forcing term. In particular, the proof 
requires a delicate analysis of the problem as well as gradient estimates for the 
fundamental solutions which can only be represented via so called Fox iL-fucntions. 
In the analysis of these special functions we use number theoretic tools to obtain 
their behavior up to the first derivatives. A particular difficulty in all the analysis 
is caused by the fact that the fundamental solutions Z and Y have singularities 
also for positive times. This causes integrability problems and requires a delicate 
analysis. 

We continue to study decay results by two additional approaches. In the first one, 
we use Fourier techniques to build optimal L^-decay estimates for mild solutions 
of the aforementioned Cauchy problem. Contrary to the standard caloric func¬ 
tions, the decay rate does not improve with high enough dimensions, but there 
exists a critical dimension at which the decay rate of bounded domains is achieved. 
This critical dimension phenomena is brought by the introduction of the fractional 
Riemann-Liouville time-derivative and such behavior is not observed in the case 


REPRESENTATION AND ASYMPTOTIC BEHAVIOR OF SOLUTIONS 


3 


a = 1. This also substantially complicates the analysis and we are required to use 
Riesz potential estimates to obtain the decay results. Thus the theory is markedly 
different from that of the standard heat equation. 

Finally, we turn into studying the decay of weak solutions where we can consider 
operators C with general measurable kernels. We show that the L^-norm of a weak 
solution, which is defined in a variational formulation, is a subsolution to a purely 
time-fractional equation. On the other hand, the exact behavior of the solutions for 
such problems is well-known and, therefore, we may use the comparison principle 
to conclude the result - even in such a general context. While our method gives the 
optimal decay rate in the case a = 1, the energy methods used in the proof cannot 
discriminate between large and small dimensions. Consequently, we are not able 
to obtain the non-smooth decay behavior - and the consequent critical dimension 
phenomenon - with respect to the dimension. Thus, it remains an open question 
whether our decay result is optimal in this context. 

2. Preliminaries and main results 

2.1. Notations and definitions. Let us first fix some notations. We denote the 
space of /c-times continuously differentiable functions by and (7° := C. 

The Riemann-Liouville fractional integral of order a > 0 is defined for a = 0 as 
'■= I, where / denotes the identity operator, and for a > 0 as 

= 7^ / (i-'r)““V(T)dT= (g«*/)(t), (2.1) 

r(a) Jo 

where 

—1 

^ fH 

is the Riemann-Liouville kernel and * denotes the convolution in time. We denote 
the convolution in space by * and the double convolution in space and time by S. 

The Riemann-Liouville fractional derivative of order 0 < a < 1 is defined by 

af/(0 = sT-“/(i) = sfyb) - -)-"/(r)*. (2.2) 

Observe that for sufficiently smooth / and a S (0,1) 

d?{f - fio)m = =■. ( 2 . 3 ) 

the so-called Caputo fractional derivative of /. In case ex = 1, we have the standard 
time derivative. 

Let 

u{0 = = (27r)"'^/2 [ e"“''^/(x)da: 

and 

denote the Fourier and inverse Fourier transforms of u, respectively. We define the 
fractional Laplacian as 

Next we define the concept of a classical solution of (HI, with £ = (—given 
with an initial condition u{0,x) = uo{x). 
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Definition 2.4. Let 0 < a < 1 and 0 

/ G C([0, oo) X R"^). Then a function u G 
the Cauchy problem 

d^{u{t,x) - uo) + 

u(0, x) 
if 

(i) defines a continuous function of x for each t > 0, 

(ii) for every x G R'^, the fractional integral as defined in (El]), is 

continuously differentiable with respect to f > 0, and 

(iii) the function u(t, x) satisfies the integro-partial differential equation of (12.511 
for every {t,x) G (0,oo) x R'^ and the initial condition of (12. 5p for every 
X G R^*. 


< /3 < 2. Suppose Uo G C'(R'^) and 
(^([OjOo) X R'^) is a classical solution of 

= in (0,oo)xR^, 

= uo{x), in R'^, 


We remark that under appropriate regularity conditions on the data, existence and 
uniqueness of strong L^’-solutions of (1^ follows from the results in , which are 
formulated in the framework of abstract parabolic Volterra equations, see also the 
monograph [5S] . 

Next we turn in to the weak solutions to equation EU. In place of the fractional 
Laplacian we will consider a more general class of elliptic operators. In this con¬ 
text, we avoid using the Fourier transform and the corresponding definition for the 
fractional Laplacian is given by its singular integral representation 


where P.V. stands for the Cauchy principal value and c is a constant. In [33] it is 
shown that (—A)^tt(x) is a continuous function whenever u is locally in C'^(R‘^) 
and 


|u(x)| 


/R<i 


1 -I- \x\^+P 


dx < oo. 


We will study the weak formulation where we define the operator through a bilinear 
form. We begin by setting up the problem. 


We define the fractional Sobolev space IF 2 ’^(R'^) for /3 G (0,2) as 


:= G L2(R^) : g ^2^ 

[ \x-y\ ^ 


X R'*) 


endowed with the norm 




f f [ -“(ff dxd. 


1/2 


Let 0 < A < A and define the kernel iG : R'^ x R'^ —>■ [0, oo) to be a measurable 
function such that 

|x-2/|d+/3 - - \x-y\d+0 

for almost every x, y G R'^ and for some (3 G (0, 2). Consider the bilinear form 

£{u,v):= / K{x,y)[u{x) - u{y)] ■ [v{x) - v{y)]dxdy 

JR'i jR'i 


(2.7) 
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for any u,v G Let ip € We now define an elliptic operator C 

by 

{Cu,(f) = £{u,(p), 


where (•, •) stands for the duality pairing on V' x V with V = Observe 

that if K{x,y) = c{d, P)\x — the operator £ defined here gives the fractional 

Laplacian of (EH). 

We study the Cauchy problem for weak solutions of the equation 

df{u — uo) + Cu = 0. (2.8) 

In the case a = 1, a weak solution is defined in the classical way. Letting T > 0, 
a natural parabolic function space for defining weak solutions on [0, T] x in the 
case a G (0,1) is given by 

F„{T) := {v G such that 

gi_a * u G C'([0,T];L^(M‘^)) and (5i_„ * i;)|t=o = 0}, 

cf. gZ]. 

The definition of weak solution (in the case a G (0,1)) is now the following. 

Definition 2.9. Let uq G L^(R'^) and u : [0,oo) x ^ R be such that for any 
T > 0 we have u|[o,t]xR‘^ £ Fa{T). Then we say that u is a weak solution of 
equation (12.81) with initial condition rt|t=o = uq if for all T > 0 

— [gi-oL * iu{t, x) — uo{x))]dt(p{t, x) dx dt 


10 Jr’>' 


K{x,y)[u{t,x) - u{t,y)] ■ [(fit,x) - ipit,y)] dxdydt = 0 


Jo Jr'^ Jr'' 

for all test functions ^ G W^’^{[0,T]] L'^{R‘^))r\L^{[0,T]-,W^''^ 
in 


with p\t=T = 0 


We recall that existence and uniqueness of weak solutions in Fa{T) has been studied 
in |47j . even in a more general context. 


2.2. Fox iJ-functions. The Fox iJ-functions are special functions of a very gen¬ 
eral nature and there is a natural connection to the fractional calculus, since the 
fundamental solutions of the Cauchy problem can be represented in terms of them. 
Since the asymptotic behavior of the Fox iJ-functions can be found from the liter¬ 
ature, the Fox iJ-functions have a crucial role also in our asymptotic analysis. We 
collect here some basic facts on these functions. 


Let us start with the definition. To simplify the notation we introduce 


.— ((ufc, Ofc), (Ufc^i, ), . . . , (Up, tTp)) 

for the set of parameters appearing in the definition of Fox iJ-functions. The Fox 
iJ-function is defined via a Mellin-Barnes type integral as 


where 


Tjmn 

^pq 


iz) 


TTmn\ \ 1 

J 


1 

27rz 






njLi b(&j + Idjs) lli=i r(l -aj- a^s) 
riLn+i r(ai -b ais) U%m+i r(l - bj - Pjs) 


( 2 . 10 ) 


( 2 . 11 ) 
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is the Mellin transform of the Fox H-function iF™" and £ is the infinite contour in 
the complex plane which separates the poles 

bji 


b,i = -^—- (j = / = 0,1,2,...) 


of the Gamma function r{bj + Pjs) to the left of £ and the poles 


1 - a,; 


^ik — 




(z = 1,.. .,n; fc = 0,1,2,...) 


to the right of £. 

We will need the following properties from Chapter 2 of [2U] . 
Lemma 2.14. Properties of Fox Fl-functions: 

A) jL frl 1 _ —r I (0.1), (ai,ai)i,p] 

(V dz^pg L^l (bj,/3j)i,, J - ^ "p+1.9+1 (bj./3j)i,,, (1.1) J 

Aq ) Ijmn r_ —11 1 rpTim r_| (l“^j 1 

m J^pq F I (b,-,/3,)i,, J - "gp FI (l-ai.ai)i,pj 

A^.j) i.p — 1 . — r I (o^ .CKi) i.p_i 1 

(lllj tlpq J - -”p-l,g-lFI (b,-.ft)2,, J- 

Aq.) OCT rrmn Fj-p I (“t .“i)l,p 1 p— ct rpCTi.’T'+l Dpi (O.p), (ai,ai)l,p] 

F I J - .«p+l,g+lLi: I (ct.p) j 


(v) For b > 0 and x > 0 we have 

nCO 

— — \h 9 '^ T~'^ \ (“i.aili.p. (1-^+7.^)l 

- ^ "p+2.g 2 ; I J 

All) r PJ™'" lyl Fi.“i)i,p 1 _ prmra f I (CTi+“i,“i)i.p 1 

(my zUpq L^l J -Pipg FI (b,+ft,/3,)i,, J- 


( 2 . 12 ) 


(2.13) 


Proof. The first four properties are straightforward calculations based on the Mellin- 
Barnes integral representation of Fox iJ-functions. For properties (v) and (vi), we 
refer to Corollary 2.5.1 and Property 2.5 of [53], respectively. □ 

Remark 2.15. There are some restrictive conditions on the parameters appearing 
in (u) (for details, see [531 Corollary 2.5.1]). The conditions are required for the 
convergence of the integrals. Since (u) represents a Hankel transform formula for the 
Fox H-functions and the Fourier transform can be written as a Hankel transform, we 
will use (u) in the proof of Theorem 2.22 to calculate the inverse Fourier transform 
of ^ I—t |'?|^^(£^), which is not integrable in general. But since both sides of (u) 
depend analytically on our choice of parameters, the identity 

(/,£) = (/,£), (2.16) 

where (•,•) denotes the duality pairing on <S' x 5 with S denoting the space of 
Schwartz functions on allows us by analytic continuation to conclude that 
(v) is valid even for a wider range of parameters than the range given by the 
restrictions in [^ Corollary 2.5.1]. Then the Hankel transform formula (v) has to be 
understood as the generalized Fourier transform (12.161) . For details on generalizing 
integral identities we refer to [23] . 
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An important special case of the function with the parameters (a^, ai)ip = 

(0,1) and {bj,l3j)i^2 = ((0,1), (1 — a,/?)) is the two-parameter Mittag-Leffler func¬ 
tion 

CX) 

i?a./3W = gr(/3 + afc)- 

It appears in the fundamental solutions of the Cauchy problem for integro-ordinary 
differential equations. Since the problem (12.51) formally transforms into 


-uoiO) 


u{0,O = uo{^), 


(2.18) 


the fundamental solutions in the Fourier domain can be formally expressed in terms 
of Mittag-Leffler functions. It can be shown rigorously that in our case the funda¬ 
mental solutions Z and Y satisfy 

(2.19) 

and 

Yit,0 = (27r)-''/V-i£;„,„(-|C|''t“). (2.20) 

The Mittag-Leffler function Ea^a{—x) is known to be completely monotone for 
X € K.+ and it has the asymptotics 

Ea,a{-x) ^ 2 ’ ^ ^ ®+- (2-21) 

1 + 

For Aq,,! we have the asymptotic behavior 

Ea^i{-x) ^ a;eR+. (2.22) 

1 -I- a; 


The asymptotic behavior (I2.21|) follows from an integral representation 


Ea,^ (-^) 


1 f , 

-— / - dt 

2tti Jc t°‘ — z 


where C is an infinite contour in the complex plane. For details we refer to EH 
Chapter 18]. Alternatively, one can use the connection to the Fox iL-function and 
use the asymptotic behavior known for the Fox H functions, see Section [3^ For 
the function Y we obtain the asymptotics 


yit.o 


^a-l 


(2.23) 


2.3. Main results. Our first theorem states that a classical solution of (E^ has 
an integral representation involving Z and Y mentioned above. Since we de¬ 
fined the fractional Laplacian via the Fourier transform, we need to guarantee 
that .^)) determines a continuous function in x. In particular, by the 

Riemann-Lebesgue Lemma, this is true if | • \^u{t, ■) G L^(R‘^). We also need that 
u{-,x) is a continuous function up to 0 for all x G For these purposes, we 
impose the condition 

\f{t,0\<Cm)\ (2.24) 

for the forcing term /, where the function g satisfies 

(l + |.|/3)g(.)GLi(R'^), 


and C > 0 is a constant which is uniform in time. 


(2.25) 
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Theorem 2.26. Letuo G be such a function thatuQ G and let f be a 

function satisfying f{t, •) G for all t >0 and (I2.24|) with g satisfying (12.251) . 

Define 

(2.27) 


Z{t,x) = 'K iJgl[2^t“|a; 


|-/3| (0.1), (o.-l)- 

I (0.1), (O.a) 


and 


Y{t,x) = 7r-‘^/h‘^-^\x\-‘^Hli[2^nx\-P\ (2.28) 
Then the function 

'l’{t,x)= Z{t,x-y)uo{y)dy+ / Y{t - s,x - y)f{s,y) dy ds 
jR‘i Jo JR'^ 

is a classical solution to problem (ESI). 

Remark 2.29. In our asymptotic analysis we prefer to use the similarity variable 
R = t~°‘\x\^ similarly as in [5D]. Therefore, it is desirable to use the property {ii) 
of Lemma 12.141 and write Z in a form 

Z{t,x) = TT-'^/'^\x\-’^Hll[2-H-°‘\xf\ (I’l)] (i,/ 3 / 2 )] (2-30) 

and y in a form 

Y(t,x) = (i|(;/ 2 ), (i,’i“)! ( 1 ./ 3 / 2 )]- (2-31) 


Observe that in the special case /3 = 2, we obtain the time-fractional diffusion 
equation. Its decay properties have been studied in and for the behavior of its 
fundamental solution, we refer to [32] . If we restrict our formula (12.271) to the case 
/3 = 2, it reduces to 




(2.32) 


Using the properties (ii) and (hi) of the Fox i7-function from Lemma [2.141 gives 


H^l[4t° 


-21 ( 1 -#, 1 ). ( 0 . 1 ). ( 0 . 1 )] _ 


(0.1). (O.a) 


= [4t“r 


-2| (1-fl). (0.1)- 

(O.a) 

2,-a| (l.“) 


- ff20r^m|2.-a| lUaJ ] 

- ^12 141-^1 I (|,1), (1,1) J- 


Therefore the formula (12.3211 reads as 

Z{t,x) = TT-'^/^\x\~‘^H^^[^\x\'^t 
which is exactly the same as obtained by Kochubei in [32l Formula (18)]. 


2 + -a| (l.«) ] 

(i.l). (1.1) J’ 


As explained earlier, the functions Z and Y can be derived by taking the Fourier 
transform with respect to the spatial variable x and the Laplace transform with 
respect to time in dLH). For more details we refer to m- Our contribution is in 
showing that they induce a representation formula, even for relatively rough initial 
and forcing data. 


Adopting the notion of the Green matrix from [20], we call the pair {Z,Y) the 
matrix of fundamental solutions of equation (ESI). Next we define the concept of 
mild solutions by means of the above representation formula. 


Definition 2.33. Let ug and / be Lebesgue measurable functions on and 
[0, 00 ) X R"^, respectively. The function u defined by 


i{t,x)= / Z{t,x-y)ug{y)dy+ / / 

JR ‘1 Jo JR-^ 

— ■ '^initit^x') -\~ UpQrc{t^X^ 


- s,x- y)f{s,y)dyds 
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is called the mild solution of the Cauchy problem (12.51) whenever the integrals in 
the above formula are well defined. 


We are particularly interested in the case where the data belong to some Lebesgue 
spaces. Note that our case differs from the usual heat equation. For example, in 
the case of the heat equation it is enough that uq € C'(R'^) n for the above 

defined u to be the classical solution of the homogeneous equation. As we shall see, 
for d > 2 and a < 1 the function Z{t, x) has a singularity not only in t, but also in 
X, which implies that more smoothness on mq is required. The function Y also has 
a singularity both in t and x. Notice that this resembles the Laplace equation, for 
which the fundamental solution u{x) = c{d)\x\'^~‘^ has a singularity at a; = 0. In a 
sense this reflects the elliptic nature of the nonlocal PDF when a < 1. 


Next we turn in to the decay of mild solutions. We give a quantitative rate at 
which the solution decays to its fundamental solution and, moreover, if the first 
moment of the initial datum is finite, we can say even more. These results are 
analogous with the ones for the heat equation in m- However, unlike in the case 
of caloric functions, we need to restrict our study of the L^’-decay to a certain range 
of possible values of p. This is caused by the fact that the fundamental solution 
lacks integrability for large enough p. Note that this does not happen for the heat 
kernel, which belongs to for all t > 0. In the limiting case we prove a 

convergence result in the weak L^-norm. 


Denote 


and 


Ki(/3,d) 


d 

d-/3+l’ 
oo, for 


for d > /? — I, 
d < /? — I 


K2(/3, d) 


d-2/3’ 


if d > 2/3, 


otherwise. 


In order to obtain decay for the solution, we need to assume that there exists a 
7 > 1 such that 

ll/(i,-)llLqR<^)<(l + 0“^ ^>0. (2.34) 


Set also 




uo{y)dy and Mforc 




f{t,y) dy dt. 


With this notation we have the following result. 

Theorem 2.35. Let d > 1, uq G L^{W’‘) and f S x Suppose f satis¬ 

fies (12.341) with some 7 > 1. Assume that u is the mild solution of equation (12.51) . 

(i) Then 

t'f ^^~p')\\uinit{t, •) - MinitZ{t, •)||iP -)> 0 , flst -)• 00, 
for all p € [ 1 , Ki), and 

t^^^^^~^'^~°‘\\uforc{t, ■) - MforcYit, OIIlp 0, t ^ OO, 
for all 1 < P oo, if a = 1 or d < 2/3, and for p G [1, K 2 ), if d > 2/3. 

(li) Assume in addition that |||a;|Mo||Li < 00 . Then 

t~ ^^~p')\\Uinit{t, ■) - MinitZ{t, ■)||lp < t~'^ , t > 0. 

Moreover, in the limit case p = d) we have 

q (/ 3 - 1 ) _a 

i ^ (^; *) MijiiiZ{t, •') II 00 t ^ , 


t > 0 . 
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Continuing on decay results, we now turn to study the L^-decay of mild solutions. 
Observe the critical dimension phenomenon that the decay rate does not improve 
when the dimension is increased after d > 2/3. Thus, the non-local case is markedly 
different from that of the standard caloric functions. Importantly, in Section [7] we 
will also show the decay rate provided here is optimal. In particular, the decay rate 
below is sharp for all initial data uq such that uq dx ^ 0. 

Theorem 2.36. Let a. G (0,1), d > 1 and d ^ 2/3. Suppose u is the mild solution 
of the Cauchy problem (12.51) with uq € L^{W^) n and / = 0. Then 

I|w(^,-)IIl= i>0. 

Moreover, in case (3 = 2/3 we have 

\\u{t,-)\\L^,^ <t-°^, t>0. 


Finally, in the following theorem we turn in to the decay of weak solutions. The 
proof is based on a comparison principle and a priori estimates. It is an open 
question whether the decay rate here is optimal as it is not as good as the one 
obtained by the Fourier methods in the previous theorem. The same phenomenon 
is present already in the case of the time fractional diffusion . Observe that our 
method gives the correct decay when applied to the heat equation. 

For s € (0,1) we set 






\v{x) - v{y)\ 
\x — y\d-+^ 


dx dy, 


which is the Gagliardo-seminorm of the Sobolev Slobodecki space 


Theorem 2.37. Let uq G and suppose the kernel K satisfies dZIl) 

with some /3 € (0,2). Let u be the weak solution of equation (12.811 with initial 
condition u\t=o = uq, and assume that 


Then 


f \u(t, •)! /3 dt < oo for all T > 0. 

Jq ^V2A(Rd) J 


(2.38) 


■“(^i < (1 +1) , t > 0. 


Remark 2.39. (i) As our proof shows. Theorem 12.371 ('trivially') extends to the case 
where the kernel K also depends on time t, that is K = K{t,x,y), provided that 
K is measurable on (0, oo) x x and (12.71) holds a.e. in this set with K{t, x, y) 
in place of K{x, y). In this more general formulation, our result can be also applied 
to certain quasilinear equations which satisfy suitable structure conditions. 

(ii) The authors believe that by careful estimates for appropriate approximating 
equations (as in [47]) combined with Gagliardo-Nirenberg inequalities one can show 
that the weak solution of equation (12.81) always satisfies the technical condition 
(12.381) provided that uq € T^(R‘^) nL^(R'^). For the sake of simplicity we do not go 
into the details here. 


3. Auxiliary tools 


We recall some classical results which are needed in the theory. 
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3.1. Review of harmonic analysis. Let f *g denote the convolution of f,g on 
We recall the Young’s inequality for convolutions: for any triple 1 < p,g,r < oo 
satisfying 1 + i = i + i 

||/*5l|L^<||/||Lr||ff||L., / e e (3.1) 

We also recall the strengthened version for weak type spaces: Let 1 < p, g, r < oo 
satisfy 1 + ^ = ^ + i. Then 

j r p q 

\\f*g\\Lr<Cip,q,r)\\f\U.,^\\gU,, / e g G (3.2) 

see [231 Theorem 1.4.24]. In the case g = 1 there also holds 

\\f*g\\LP.^ < C{p)\\fh..^\\gh^, f G L^’(R‘'), g G ^^(R'^), (3.3) 

for all 1 < p < oo, see [Ml Theorem 1.2.13]. 

For the nonhomogeneous problem we need the integral form of the Minkowsky 
inequality in the following form. Let 1 < p < oo and F’ be a measurable function 
on the product space R+ X R'’*. Then 

if if \P{tT^)\dx\ dt]” < [ ( f \F{t,x)\^ dt]” dx. 

We will also need the following decomposition lemma from m- 

Lemma 3.4. Suppose f G L^(R^) such that \x\\f{x)\dx < oo. Then there 
exists F G L^(R‘^;R‘^) such that 

f= ^o + divF’ 

in the distributional sense and 

ll^llLi(R<i;R<i) <Cd f |a;| \fix)\dx. 


We will also need the boundedness of the Riesz potential 


(-A)-^/:= cd,^ 


/(y) 


Vd \x - y\^-d 


dy, 


for 0 < /3 < d. We have the Hardy-Littlewood-Sobolev theorem on fractional inte¬ 
gration |241 Theorem 6.1.3]: 

Theorem 3.5. Let 1 < p < d/fi and f G LP(R^). Then 

ll(-A) ^ ^ ll/llLP(R'i) 


for p > 1 and 


dp 

d — pfF 


ll(-A)-^/|| 


Ld-Pfi ’°“(R<i) 




In case p = 1, we have 
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3.2. Asymptotic behavior of the Fox iJ-functions. When developing the as¬ 
ymptotic behavior of the fundamental solution, we need the following representation 
formulas for the Fox iJ-function Here we have omitted the parameters of the 
Fox iJ-function and i?!! refers to the Fox i7-function appearing either in 12.301 or 
in 12.311 The following results hold for both functions. 

Theorem 3.6. Let either P > a and z ^ Q, or a = P and 0 < |z| < d with 
S = Then the Fox H-function is an analytic function of 

z and 

2 oo 

^23 (^) = E E , (3.7) 

i=i1=0 

where bji are given in (j2.12ll . 

The asymptotic behavior of i?!! ('^) as z —?> 0 follows immediately from (13.71) in the 
case /3 > a by calculating the residues. If 0 < a < /3 and |argz| < 7r(l — ^), then 

2 oo 

'--EE Res.=-b,, [nl\{-s)z^] , (3.8) 

i=i 1=0 

when z —> 0. Again, the asymptotic behavior follows immediately by calculating 
the residues. 

The asymptotic behavior at infinity is more complicated to derive. For details we 
refer to [5] and [121 Sections 1.3 and 1.5]. 

Theorem 3.9. The asymptotic expansion at infinity of the Fox H-function Hf^^z) 
has the form 

OO 

Hll{z) (3.10) 

fc—0 

where the constants hk have the form 
hk= lim [-{s - aik)'H 2 l{s)] 

s^aik 

(_l)fc r(6i + (1 - ai + fc)^)r(b 2 + (1 - ai + k)^) (3-11) 

" k'.ai r(a2 + (l-ai+A:)^)r(l- 63 -(l-ai + A:)^) 
in view of the relation 

^ess=a,,[nll{s)z-^] = hfeZ-“^'= = 

4. Behavior of the fundamental solutions 

We start by showing some basic properties of the fundamental solutions Z and Y. 
The first lemma provides an important connection between the functions Z and Y. 
Note, in particular, that Z and Y are identical in the case a = 1. 

Lemma 4.1. The fundamental solutions Z andY of equation (El) are connected 
via Y = dl-°‘Z. 

Proof. Observe hrst that 

9i-“/(at) = ai-“(cli-“/)(at) 

for a sufficiently smooth function / and for a constant a £ 1R+. 
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Now we combine this with Lemma [2.141 ('iv') to obtain 

Ql-arrl2r<^/3.a, |-/3| (l-f,/3/2). (0.1), (0,/3/2) n 
Ot t |x| I j 

_ .a-lrrl3 1 - 9 / 3 .ai 1-/31 (0,a). (l-f./3/2). (0.1), (0,/3/2) n 

-t t |a;| | (i-a.a) J’ 

We need to study the Mellin transform of the above Fox iL-function. That is 

r(s)r(i - as)r - f s) r(i - s) 


Hll(s) = 


r r(l — as)r(a — as) 

r(s)r(f-fs)r(i-s) 


(f s) r(Q; — as) 


= ^32(5) 


We obtain 

Y{t,x) = n-'^/^\x\-H‘^-^Hll[2Pnx\-^\ 
= dl-<^Z{t,x), 

as required. 


□ 


Before moving into providing the exact behavior of the fundamental solutions Z 
and Y, we give the following remark. 

Remark 4.2. Observe that the functions Z and Y are both non-negative and, more¬ 
over, Z induces a probability measure. 

Indeed, by Bochner’s Theorem the non-negativity follows from showing that the 
Fourier transforms Z{t, •) and Y{t, ■) are positive definite on [3]. Recalling that 
Z{t, •) and Y{t, •) can be represented in terms of the Mittag-Leffler functions E’a.i 
and for the positive definiteness it is enough to show that the functions 

/(r) = Ea^i{—t°‘r^) and g{r) = Ea^ai—t°‘r^) are completely monotone on R+ [38l 
Theorem 3]. But since the functions x i—)• Ea^i{—x), x i—>• Ea^a{—x) and x i—>• cx^~^ 
with c > 0 and f3 < 2 are known to be completely monotone on M+ [35] , we obtain 
the result. 

Finally, by (12.191) we have 

/ Z{t,x)dx = Z{t,0) = Ea^i{0) = 1, 

for every t > 0, which yields that Z{t, •) > 0 induces a probability measure on W^. 


When proving the decay estimates we will need the following asymptotic estimates 
for the fundamental solutions. We begin by studying the function Z. 

Lemma 4.3. Let d S 0 < a < 1 and 0 < /3 < 2. Denote R := |x|^t““. Then 
the function Z has the following asymptotic behavior: 

(i) If R < 1, then 

{ ^-ad/p^ Q, _ or 13 > d and 0 < a < 1, 

t““(| log(|x|^t““)| -f 1), if (3 = d and 0 < a < 1, 

if0</3<d and 0 < a < 1. 


Z{t, x) 
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(ii) If R> 1, then 

if 13 <2. 

In the special case f3 = 2 there holds 

Z{t,x)<r\x\-‘^-^. 

Proof, {i) R < 1: We start with the case 0 < a < 1. Since the asymptotic behavior 
depends on whether f3 > a or (3 < a, we have study different subcases. First of all, 
recall the definition of Z as 

z(f,x)=„ j,]. 

In order to figure out the asymptotic behavior of Z, we need to study the above Fox 
Fl-function. As it was mentioned in Section 13.21 the asymptotic behavior follows 
by calculating the residues. We provide the details for the reader’s convenience. 

The subcase f3 > a: We have 

2 oo 

= Ress=b^i (4-4) 

i=i /=o 

by Theorem 13.61 Recall the definition of the Mellin transform 

Rjlis) ■.=M{Hll[z\ (i.|)])(s) 
r(f+ |s)F(i + s)r(-s) 

F(l + as)F(-fs) 

In light of (14.411 . the asymptotic behavior is determined by the largest value of s, 
which is a pole of TL‘^{s). Now, for 0 < a < 1 the above Mellin transform has poles 
at s = —1 and s = —d/f3. Suppose first that f3 > d. The only value for d to this 
happen is d = 1, when 1 < /3 < 2, whereas for 0 < ,5 < 1 there is no such d. Then 
the asymptotics is determined by the pole at — I/^S and the behavior of the 
near zero is ~ This yields 

Z{t,x) ^ 

as required. 

Assume next that a < f3 < d. Then the largest value of s such that the Mellin 
transform has a pole is s = — 1 and we obtain Hf^{z) ^ z. This produces 

In the case (3 = d the Mellin transform has a second order pole at s = —1. Then 
the residue can be calculated as 

Ress=_i [^ 23 ( 5 ) 2 “'’] = ^hm^ ^[(s + 

= z lim -^[{1 + sfnllis)] 
s-^-l as 

+ \im U1 +senilis)-^{z-^)]. 

Since (I + s)F(I + s) = r(2 + s) and |(1 + s)r(| + ^s) = F(| +1 + §s) are analytic 
at s = —1, the limits 

lim [(1 + 3 ) 2 ^ 23 ( 5 )] and lim (1 + 
s —»■—1 as s—S -—1 
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exist. Moreover, since 

lim = zlogz, 

s —>■—1 as 

we may conclude that in this case H 2 ^{z) ^ zlogz and thus 

Z{t,x) ~t““(|log(|a;|t"“)| + l), 

again as required. 

The subcase f3 < a: Since we are interested on the asymptotics of the Fox H- 
function for z S K+, the asymptotics is given by (13.81) . Because 0</3<Q!<l<fi, 
we have d/j3 > 1 and the leading term is determined by 

Therefore 

Z{t,x) ^ 


In the special case a = 1 we see that the Mellin transform of H^Hz) reduces to 

.,21. . r(f + |s)r(-s) 

H23(^) - • 


Therefore the asymptotics is given by the pole ats = —^ = —Proceeding as 
above we end up to the desired estimate. 


in) R > 1'. We use the asymptotic behavior of the Fox i/-functions provided by 
Theorem EH 

OO 

fc =0 

for constants hk defined in (13.1111 . We aim to find the smallest value of k such that 
hk ^ 0. Let’s first study the case 0 < /3 < 2. Now 


r(i)r(i) 

F(i)r(o) 


= 0 and hi 


F(f + |)F(2) 

r(i + a)r(-§)^ 


(4.5) 


Therefore the leading term in the expansion (I3.10|) is hiz ^ so 


^ ^ ^ ^ OO 


and we obtain the claim of the lemma. 


If /3 = 2, we see from (SH) that ho = 0 and h(l) = 0, since 

r(-f) = r(-i) = oo. 

Therefore the claim is true also in this case. However, we can continue to deduce 
hk = 0 for all k G Z+. One can prove that now actually Z{t^x) decays in terms of 
R — >■ OO, but we do not need that fact in our considerations. □ 


The next lemma gives the behavior of the fundamental solution Y. The proof is 
similar to the previous lemma. 

Lemma 4.6. Let d>l, 0<Q!<1 and 0 < /3 < 2. Denote R := \x\^t~°‘. Then 
the function Y has the following asymptotic behavior: 
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(i) If R <1, then 

, if d> 2fi and 0 < a < 1, 

log(2“^|x|^i““)|, if d=2(3 and 0 < a < 1, 

^ Q, = or d < 2/3 and 0 < a < 1. 

(ii) If R> 1, then 

Y{t,x) if0<l3<2. 

In the special case (3 = 2 there holds 

Y{t,x)<t^°‘-^\x\-'^-'^. 

Proof. The proof is similar to that of Lemma 14.31 We omit the details. Once again, 
notice that in the special case /3 = 2 the function Y has indeed exponential decay 
as i? —)• oo but we do not need that fact in our calculations. □ 


Next we turn to study the behavior of the derivatives of Z and Y. 

Lemma 4.7. Let d > 1, 0 < a < 1 and 0 < /3 < 2. Denote R := \x\^t~°‘. Then the 
derivatives of the fundamental solution pair (Z,Y) have the following asymptotic 
behavior: 


(i) For the function Z we have 

\S/Z{t,x)\ - if R<1 

|VZ(t,x )|if R>1. 


(i) For Y we have for R < 1 that 


\VY{t,x)\ 


and 


t “ ^|a;| 1+2/3^ if d+ 2 > 2/3 and 0 < a < 1, 

< t“““^|a;log(2“^|a;|^t““)|, if d + 2 = 2/3 and 0 < a < 1, 

T cx(d+2) 

t°‘ 3 |x|, if a = 1 or d+ 2 <2/3, 


|vy(t,x)| if R>1. 


(ii) In addition, for the time derivative of Y we have 
\dtY{t,x)\t~^Y{t,x), if i? < 1, 
\dtY{t,x)\^t^<^-^\x\-'^-^, if R>1. 


Proof. We provide the calculations only for the gradient of the function Z. The 
other cases are handled similarly, but we omit the details. 


Recall the expression for the fundamental solution Z: 

First of all, we use Lemma [2.141 ('ii') to write the above Fox id-function as 


HM[2Pnx\-P \= id|i[2-/^t-“|j 


(14). (l.a) 

(f.f). (1,1). (1,1). 


According to Lemma 12.141 (i), we have 


_d 

dz 



(1.1). (i.«) 

(|.f). (1.1). (1.1) 


_ ^-lH-22r.,| (0.1). (1.1). (l.a) 

-Z ^^34 L^l (I £), (1,1), (1,1), (1,1) 
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Using the product rule for differentiation, we may now calculate 

For simplicity, we have here omitted the set of parameters inside the Fox H- 
functions. Next we analyse the above Fox iJ-functions by studying the correspond¬ 
ing Mellin transforms. We have 


_ ^ r(f+ |s)r(i + s)F(i-s) _ ^ r(f + fs)r(i + s)r(-s) 


F(l-has)r(-|s) 


F(l-has)r(-|s) 


= _2 r(f+ |s )r(i + s)n-s) 


2 2 


F(l + as)r(-§s) 


= -2 


r(rf^ + |s)r(i + s)F(-s) 


r(i-has)r(-|s) 

— ^^23 1*1 (^+2 


(1,1), (1,1) J- 

Thus we obtain 

'’-Z(t.z) = [2-'>(->|x|->|, W'.,, ] 


dxi 


and 


|VZ(t,x)| = 27r-''/2|ai|-‘'-i Hll[2-h-^\x\^\] 


Now the result follows from the behavior of the Fox Ff-functions. 


□ 


5. Representation formula for solutions 


5.1. Proof of Theorem 12.261 We are now ready to prove Theorem 12.261 which 
justifies calling (Z, Y) the matrix of fundamental solutions for the equation (12. 5|) . 


Proof. We need to show that the function 
'l’{t,x)= Z{t,x-y)uo{y)dy + 

=: 'i>i{t,x) + 'ii 2 {t,x) 

is a classical solution to equation dm. We divide the proof into three steps as 
there are three requirements in the definition of the classical solution. 


0 JRrf 


Y{t- s,x- y)/(s, y) dy ds 


Step I: First we need to prove that ^)) ^ continuous function with 

respect to x for each t > 0. The representation (12.1911 and the asymptotic behavior 
of the Mittag-Leffler function given by (I2.22|) give 


\z{t,0\< 


c 

1 -b 


for f > 0. Thus \£,\dZ{t,^) is bounded for all f > 0 and by using the assumption 
that Mo G we obtain 




In order to estimate F, we use the assumption (12.2411 to obtain 


\f(.s.o\<m)\ 


(5.1) 
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with S L^{W^). Combining this with (12.231) yields 

*f){t,0 = lei"" r y{t --ms, •) ds 

Jo 

l+*'(t_5)2a|^|2/3 

which establishes that | • $ f)(t, •) S L^. This together with (15.11) gives, again 

by the Riemann-Lebesgue lemma, that 

is a continuous function, as required. 


Step II: We proceed as in Section 5.3 of [20]. By Lemma [2.141 (iv) we have 


dtHMit 


12 ua| (1-1.1). (0,1). (O.f)i 


(0,1), (O.Q) 




(O.c), (l-#,f). (0,1), (O.f) 
(0,1), (0,a), (q:,Q:) 


A detailed study of the asymptotics similarly as in Lemma 14.31 can be used to 
show that d^Z{t,-) is integrable for all t > 0. Thus for all x G the function 
is continuously differentiable with respect to time for all t > 0. We 
now turn to study 'I' 2 . 

Let V := Observe that, after changing the order of integration and a 

change of variables. Lemma [4.11 gives 


v{t,x)=[ j [ Y{t - X,x-y)f{X,y)dydXdT 

Jo r(l — a) Jq J^d 

^ f f rn 1 f Y0'-X,x-y)f{X,y)dydTdX 

Jo Jx “ Cl) JR'* 

= f f Z{t-T,x-y)f{T,y)dydT. 

Jo Jr‘‘ 

Using Remark 1321 to deduce that Z is a probability density gives 


-[u(t + h,x) - v{t,x)] 




t-\-h 


-J 


/ Z{t + h - s,x - y)f{s,y)dsdy 


“i 


^ [Z{t + h - s,x - y) - Z{t - s,x - y)]f{s, y) dy ds 
0 Jr^ ^ 

t-\-h 


pt-\-h 


p pt+n 

/ Z{t + h-s,x-y)[f{s,y)-f{s,x)]dyds + f /(s,x) 
Js.'^ Jt 


ds 


- [Z(t + h - s,x - y) - Z{t - s,x - y)]f{s, y) dy ds. 
Jo Jk.'^' “ 


At this point we use the conditions (12.241) and (12.251) to give some regularity for 
the right hand side /. Our conditions guarantee that f{t, ■) is a Holder continuous 
function with the Holder exponent less than min{l,/3} uniformly in t, i.e. we have 
an estimate 


\f{t,y)- f{t,x)\<C\x-yp, t>0, 
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for any 0 < 7 < min{l, /3}. Using this we obtain 

Jt 


< 


/ Z{t + h-s,x- y)[f{s,y) - f{s,x)] dyds 

f / Z{s,x - y)\f{t + h - s,y) - f{t + h - s,x)\dyds 
Jo JR‘i 

/ f Z{s,x - y) ■ \x - yl'^ dyds 

Jo Ji 

■f 

Jo 


'0 J {\x — y\>s°^/ 

/ Z{s,x-y)-\x-ypdyds. 

lo J {\x—y\<s°‘^^} 

and continue by using Lemma 14.31 to conclude that 


nh n nh pOO 

f / Z{s,x — y) ■ \x — yP dy ds <-h s“ / dr ds 


0 

h 

< i ds < 0 


J 

Jo 


as > 0. Utilizing Lemma [UHl similarly in the second integral gives 


f 


Z{s,x — y) ■ \x — j/P dy ds < —>■ 0 

/O J {\x—y\<s°‘/^} 

as h ^ 0. Here one needs to check different cases depending on the values of a, /3 
and d. 


Altogether we have that 

v{t + h,x) — v(t, x) 


lim 

h^O 


h 


= + 


0 Js.'^ 


dZ{t — s,x — y) 


f{s,y) dy ds 


and, therefore, the function “dt is continuously differentiable with respect to t. 


Step III: We need to prove that the function ih satisfies the integro-partial differen¬ 
tial equation. Our assumptions on / and the asymptotic behavior of Y guarantee 
that Y{t — ■,x — •)/(■) ■) € A^((0,t) X Therefore 

x) ^ 0 , as f — 0 , 

which means that 

d^'Ii2{t,x) = d^{'Ii2{t,x) - ^'2(0, a;)). 

Notice that as a by-product of Step II we obtained 

d°'i> 2 {t,x) = f{t,x) + f f — ylff^s^y^dyds. (5.2) 

Jo JRd ot 

We will show that (9“ -I- (—A)^/^) 4 i 2 (t, a:) = f{t,x). By (15.211 it is enough to show 
that 

(_A)^/ 2 ^r 2 (t,a:) = - / f U — —fis,y)dyds. 

Jo aR<i ot 

We start by calculating 

{-Af/^Y{t,x) = 

Recall from (I2.20p that 
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Notice that in the following calculations we have to interpret the integral properly 
due to poor decay of Y{t, •) at infinity, see Remark [2.151 

Notice that y is a radial function of ^ and for radial functions we have in general 
that [24l Appendix B.5] 

nOO 

-^(/)(iei) = / /W//V^(r|C|)dr, 

Jo 

where J(d- 2 )l 2 is the modified Bessel function. For the definition, see [U]. We use 
this formula together with Lemma 12.141 (v) and (vi) to calculate 


2^ 


TT 
= TT 


o U 


(5.3) 


(14). (l.a) 

On the other hand, combining the chain rule with Lemma 12.141 ('i'l. gives 

a,Z((,x) = (»’«)]^ 

Now by studying the Mellin transform T-L^^ and using the properties of the Gamma 
function gives 

fTiaro/s^ai 1-/31 (0.1), (0.1). (o.|)n 

^43 t |ai| I (O.c). (1,1) J 

L ' ' I (1,1), (l.a) 1 

Inserting this into (15.3p yields 

dtZ{t,x) = -T-\\^fY{t,0){x). 

By the above calculation 

'■* dZ{t — s,x — y) 


T 


'0 JRd 
-t 


dt 


-f{s,y) dyds 


= f F{dtZ{t-s,0)hs,0ds = -\^f [ Y{t-s,0f{s,0ds. 

Jo Jo 

Now using the growth condition of function / (cf. (j2.24p L we have that | • Y{t — 
s, ■) f{s, ■) ds € and therefore it has a unique inverse Fourier transform. We 

obtain 

( f f Y{t-s,x-y)f(s,y)dyds 
\Jo JR"! 

' (l«tO-»,{)/(».0* 

=-/7 

Jo JR"! 


dt 


Therefore 

as claimed. 


{d^ + i-A)P/^)^SJ2it,x) = fit,x), 
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Let US now study the first integral. By using the asymptotics of Z as in iStep II, it 
is straightforward to show that 

/ Z{t,x - y)uo{y)dy ^ uo{x), as t0. 

Jr’‘ 

A similar argument as for i ['2 produces 

9t [ [ X - y)uo{y) dy - uo(a;)] + [ Z{t, x - y)uo{y) dy = 0. 

We omit the details. 

Now III satisfies the initial condition by the superposition principle. 

Step IV: Finally we have to prove that dt is a jointly continuous function in [0, oo) x 
The continuity at t = 0 is established in Step III. If t > 0, the continuity in 
both variables follows from our conditions given for uq and /, which guarantee that 
uq and / are continuous and uniformly bounded. Then the asymptotics of Z and 
Y given in Lemmas 14.31 and 14.61 together with the Lebesgue dominated convergence 
theorem imply the continuity. This finishes the proof. □ 


6. Large-time behavior of mild solutions 


We begin by calculating an L^’-decay estimate for the fundamental solution Z, 
which is given in the following lemma. 

Lemma 6.1. Let (i>l, 0 <q:<1 and 0 < /3 < 2. Then Z(t, •) S L^(R'^) for any 
t > 0 and 

\\Z{t,-)\\LHR^)<t-t(^-i\ t>0, (6.2) 

for every 1 <p < K^{P,d), where 

K3 = K3(/3, d) := I ^ ^ (6.3) 

I oo, otherwise. 

Moreover, if a = 1 or 1 = d < ft, then (16.21) holds for all p € [l,oo]. Finally, for 
d > ft and 0 < a < 1, we obtain 

||Z(t,.)|| . <t-“, t>o. 


Proof. We begin by decomposing the L^-integral of Z as 

\\Z{t, < ( Z{t,xYdx+ f Z{t,xY dx. 

J{R>i} J {R<V 

In view of Lemma El we have for all dimensions d and for all values 1 < p < oo 
that 


[ Z(t,xYdx< [ t°‘P\x\-‘^P-f^P dx 

nOO 

< / t°^Pr-‘^P-PpY-^dr<t-t(P-^\ 


and thus 


f Z{t,xY dx\ <t T p) for all I < p < oo and t > 0. 
J{R>1} / 


(6.4) 
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We come now to the estimate for the integral where R < 1. In the case a = 1 or 
P > d and 0 < a < 1, we have for all 1 < p < oo that 




Z{t, xY dx < 


t ^ dx< I t 0 ^ dr <t 0 ^ . 




2idp 


dp I ad 


i{R<n 


li B = d and 0 < a < 1 we estimate 


[ Z{t,xYdx< [ t-‘^Pi\\og{\xft-°‘)\ + lYdx 

J{R< 1 } d{R<l} 

< / f-“P(|log(/t-“)|+ 1)'’ 

Jo 

< [ t-“P+“‘'/^(|log(s'^)| + l)^s'^-^ds 
Jo 

< ^-ap+a _ ^-^(p- 1 )^ 

for all 1 < p < oo. Note that the condition P > d can only happen if d = 1. 
Finally, if 0 < /3 < d and 0 < a < 1 we have 




/{R<1} 


Z(t, xY dx 


< 


t"“P|a;r"'P+^Pda; < / dr 

' d{_R<l} Jo 

< ^-ap+|[d-(d-/3)p] ^ 


whenever the last integral is finite, that is, whenever 

P < = f^oiYd). 

d — p 

Combining the previous estimates we see that 

[ f Z{t,xYdx\ for all 1 < p < K3(/3, d) and t > 0. (6.5) 

\d{R<l} J 


Observe that by Lemma HTITI we have Z{t, •) G L°°(R) for all t > 0, provided a = 1 
or /? < d, and moreover, we have the estimate 

IIZ(t,x)llL^ 

which proves the second statement. 

For the weak-L^’-estimate we set p = We need to estimate 
\\Z{t,-)\\Lp,oo = sup^Xdz(t,x)W^ ■ A>o|, 

where 

dz{t,x)W = |{a; e : Z{t,x) > A}| 

denotes the distribution function of Z{t,x). Using again the similarity variable 
R = t~°‘\xY we have 

< 2 {\\Z{t, x)x{R<i}{t)\\LP'°° + \\Z[t, x)x{R>i}{t)\\LP'°°) ■ 

Employing (16.4L we find that 

ll^(^,ai)x{R>i}(t)||Li>.» < \\Z{t,x)x{R>i}it)\\LP < Ct~f (^"p) = Ct““. 


( 6 . 6 ) 



REPRESENTATION AND ASYMPTOTIC BEHAVIOR OF SOLUTIONS 


23 


For the term with R < 1 we use the case 0 < /3 < d of Lemma 14.31 to estimate 
dzit,x)xiR<i}{t)W = l{2; G : Z{t,x) > A and R < 1}| 

<|{a;eR'^: X < Ct-‘^\x\-‘^+^}\ 

= Kcc e R'* : |x| < {Ct-°‘X-^)^}\ 

< Cl {t-°‘x-^)^. 

This shows that 

and thus 

\\Zit,x)x{R<i}{t)\\Lp.-- < 

This finishes the proof. □ 


As a simple consequence of the above lemma we obtain the following decay result. 

Proposition 6.7. Let d>l, 0<a<l and 0 < /3 < 2. Assume that u is the mild 
solution of equation (ESI) with / = 0 and uq € L‘?(R‘^), where 1 < q < oo. Then 
the following hold: 


(i) if q = oo we have 

^ l|wo||L“(R<i), t > 0; 

(ii) if 1 < q < oo and d > q(3, we have for every r G [q, that 

•)llL'-(R‘i) ^ i > 0, (6.8) 

and if, in addition, 0<a<l<d we obtain 


\uit,-)\\ ^ 

Z, rf-9/3 ’ 




<t~ 


t > 0; 


(hi) if I < q < oo and d = qfi, the estimate (EHl) holds for every r G [q, oo); 
(iv) if d < ql3 or a = 1, the estimate (EH! holds for every r G [g,oo]. 


Proof. Let p be defined via 



r p 


1 


(6.9) 


For such p, q and r we may use Young’s inequality for convolutions to obtain 


W'^it, .)||l^ = ||^(L •) *^o(-)IIl^ < \\Z{t, ■)\\l4uo\\l.. (6.10) 

The idea is now to use Lemma | 6 T] to estimate the L^-decay of Z on the right hand 
side of the above estimate. We only need to consider different cases corresponding 
to the different choices of the parameters. 

Recall that by Lemma 16.11 we obtain 

\\Z{t,-)hHM^)<t-t(^-i'>, ( 6 . 11 ) 

for 1 < p < K 3 (/ 3 ,d), where K 3 is as in (16.31) . Now, claim (i) follows directly from 
choosing p = 1 , r = 00 and g = 00 in (16.101) . 

On the other hand, a straightforward calculation shows that for r G [g, ), we 

have 1 < p < . We again use the above L^’-estimate for Z together with (16.101) 

to obtain the claim. 
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For d> qf3 and r = j we have from (16.91) that p = . Now we may use the 

second part of Lemma [Ol to obtain that if d > 2 and 0 < a < 1, then 


\\Zit,-)\\ ^ oo <t~ 


t > 0, 

which together with Young’s inequality for weak L^-spaces gives 

||wo||L<i(R<i) 

< 


\\uit,-)\\ ^ < ||Z(t,-)|| ,o 

(R‘1) 


as required. 

For (iii), observe that inserting q = d/fi and p G [1, in (16.911 gives r G [g, oo). 
Similarly, inserting q G (d//3, oo] and p G [1, in (16.91) gives r G [g, oo]. This 
yields the first claim of (iv). If a = 1 we, in turn, by Lemma 16.11 obtain the 
LP-decay (16.111) for any 1 < p < oo and we may again use Young’s inequality, as 
in (|6.10L to obtain the claim for any r G [g, oo]. □ 


We continue by studying the above type of results for the inhomogeneous equation. 
First we need the L^-decay estimates for the fundamental solution Y. 

Lemma 6.12. Let d>l, 0<a<l and 0 < /3 < 2. Then Yit, •) G and 

lim-)llLP(R^) t>0, (6.13) 

for every 1 < p < K 2 , where 


K2 = K2(/3,d) 


d-2/3 ’ '^f d> 2/3, 

oo, otherwise. 


Atthe borderline p = K 2 , we also haveford>2P thatYlt,-) belongs to L 
and 

llY(t,.)||^^„. t>0. 

Finally, if a = 1 or d < 2/3, estimate (I6.13|l holds for all p G [1, oo]. 


Proof. The proof is similar to that of the function Z. We give the proof in the case 
d < 2/3 and 0 < a < 1 as an example. We begin by decomposing the L^-integral of 
Y as 

\\Y{t,-)\\lp = f Y{t,xfdx+ [ Y(t,x)Pdx. 

J{R>1} d{R<A 

By Lemma 14.61 we have for all dimensions d and for all values 1 < p < oo that 


'{«>!} 


Y (t, x)^ dx 


< 


t'^°‘p-p\x\-‘^p-hp dx 


I{R>1} 


^OO 

< ^2ap-p / ^-dp-/3p^d-l < ^(a-l)p-^(p-l)^ 


and thus 


[ f Y{t,xYdx 
\J{R>i} 


< ^ "3" p) for all 1 < p < OO and t > 0. (6.14) 










REPRESENTATION AND ASYMPTOTIC BEHAVIOR OF SOLUTIONS 


25 


We come now to the estimate for the integral where R < 1. Again, by Lemma 
we have 

[ Y{t,x)Pdx< [ f r‘^-^dr 

J{R<1} J{R<i} Jo 

for all 1 < p < oo, which finishes the proof of the first statement in this case. Since 
in this case even Y{t, •) G we see that the second statement holds as well. 

The weak-LP estimate is done similarly to Lemma [6.1 1 We omit the details. □ 


Again, we may use the above estimates to prove a decay result concerning the 
source term /. Here we need to impose a decay condition similar to (I2.34|) for the 
source term. We obtain the following proposition. 

Proposition 6.15. Let d>l, 0<a<l and 0 < P < 2. Assume that u is the 
mild solution of equation (12.511 with uo = 0 and f{t, •) G L'J(]R^) for each t > 0 and 
for some q G [1, oo). Assume further that f satisfies the deeay condition 

+ t>0, (6.16) 

for some 7 > 0. Then we have in case 7^1 

(i) */1 < 9 < 00 and d > g/3, we have for every r G [ 9 , 73 ^) 

IKtG)llL 7 R‘^) i>0; (6.17) 

(ii) */l < 9 < 00 and d < 9 / 3 , the estimate (16.1711 holds for every r G [ 9 , 00 ). 

In the case 7 = 1 , the assertions (i) and (ii) are valid with ( |6‘.i7| ) replaeed by 

\Ht,-)h.^^.^<t<^-^-t^-.-Jhog{l + t), t>0. (6.18) 


Proof. The proof is now an easy application of the integral form of the Minkowsky 
inequality, the Young inequality for convolutions and Lemma 16.121 

Using the Minkowsky inequality, we have 

l|u(t, OlliUR'i) < [ ([ [ Y{t-s,x-y)f{s,y)dy dx) ds 

Jo ^ 

for 1 < r < 00 . 

Similarly as in the proof of Proposition dsm, we choose p such that 

1 1 1 

1 H— — —I—. 
r p q 

Then the Young inequality for convolution yields 

ll^(^G)llL'-(R‘i) < / ll^(^ “ S; •)I|lp(R‘*)II/('S, •)IIl''(R‘^)'^®- 

Jo 

We split the integral into two parts as follows 


f \\Y{t-s,-)\\^. 

^0 

~(^J ) ll^(^ “ S; •)llLr(R‘^) ll/('S, •)llL9(R<i)<^'S =: ^1 +-^2- 

Recall that by Lemma 16.121 we have 

r(iG)llLr(M.) t>0. 


(6.19) 


( 6 . 20 ) 


( 6 . 21 ) 
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for 1 < p < K 2 where 

_ f d-2l3 ’ d > 2/3, 

1 oo, otherwise. 

If d > 2g/3, for r € [q, 3 [q, we obtain from (I6.19|) that correspond¬ 
ingly p e [1, Therefore, we may use (16.211) to estimate the L^^-norm of Y 

in (I6.20|) . On the other hand, if d < 2q/3, the different values of p G [^,K 2 ) yield 
the corresponding choices of r in [q, oo) and, thus, we may again use (16.211) to 
estimate (16.201) . 


We continue the estimate by using (16.211) . For the first integral we observe that 
I < t — s < f and hence (16.211) together with the decay condition (I6.16|l implies 


^t/2 pi 

Jo Jo 

which gives the desired estimate for Ii. 


In the second integral we need to take care of the singularity of Y{t, •) at t = 0. 
The integral converges if and only if 

This gives 1 < p < provided d > /3. If d = /3, this holds for all p G [l,oo), 
and if d < /3, this estimate is always true. Observe that this restriction gives the 
different choices of r in the items (i) and (ii) of the claim. We obtain 

pt p^l"^ 

I (l + s)"'^||F(t-s,-)||LP(Rd)ds<t“'^ / 

Jt/2 Jo 

for all 7 > 0. So I 2 decays faster than /i if 7 > 1, whereas for 7 G (0,1) we 
obtain the same decay rates. Observe also that, similarly as in Proposition 16.71 the 
restriction 1 < p < 7 ^ plays a role only if d > g/3. In this case, we obtain directly 

from (16.191) that r G [1, 73 ^)- D 


The last step towards the proof of Theorem 12.351 is the following gradient L^- 
estimate for Z. 

Lemma 6.22. Let d G Z_|_ and Ki(/3,d) be as in Section\^ Then S/Z(t,-) belongs 
to LP(W^]W^) for all t > 0 and 1 < p < Ki(/3,d), and there holds 

■)llLP(R‘i;R‘*) ^ ^ ^ *”^1 t > 0. (6.23) 

The estimate (16.231) remains valid for d = 1, /3 = 2 and p = 00 . 

Moreover, if p = Ki{f3,d), then we have that yZ{t,-) belongs to for 

all t > 0 and 

■)llLP.“(R‘i;R‘l) ^ ^ t > 0. 


Proof. The proof is very similar to that of Lemma [6.11 Let R = “ be the 

similarity variable. Let’s first divide the object of our study into two parts: 


/R'i 


\VZ{t,x)\Pdx 




\VZ{t,x)\Pdx^ 


I{R>1} 


\VZ{t,x)\Pdx=:h+l2. 
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For the first term, we may use Lemma [4.71 to get 


/i = / \VZ{t,x)\P dx< / dx 

< / ^(-d+P-i)p+<i-l 

Jo 


(6.24) 


provided the last integral is finite, that is 

1 < p < Ki{fi,d). 


For the second term we again use Lemma 14.71 and obtain 

72 =/' \S/Zlt,x)Y’dx< ( dx 

J{R>1} J{\x\>t°‘/P} 

< fOLp f ^-(d+p+l)p+d-l (6.25) 

for all 1 < p < oo. Thus we obtain the first part of the lemma. 

If /3 > d + 1, which is equivalent with (3 = 2 and d = 1, we see from Lemma [4.71 
that VZ{t, •) is indeed bounded. Therefore the second statement holds as well. 

Let now p = Ki{(3,d). Similarly as (16.6L we obtain 
||VZ(t,x)(t,-)||L.,~ 

< 2 (||VZ(t,a;)(t, •)X{R<i}(t)||Lp.“ + ||VZ(t,a;)(t, •)x{R>i}(i)||LP.~) . 

Employing estimate (I6.25P gives 

\NZ{t, x)x{R>l}{t)\\LP.°- < \\VZ{t,x)x{R>l}{t)\\LP 

< (l“p) < t~°‘. 

For the term with 7? < 1 we use Lemma 11771 to estimate as follows. 

d\vz{t,x)xiR<i}{t)\W = \{x G : \S/Z{t,x)\ > X and R < 1}| 

< |{a; G : A < Ct-°‘\x\-‘^+>^-^}\ 

= \{x G : |a;| < (C'f-“A-i) 

< Cl (t-“A-^)^™ . 

This shows that 

and thus 

\\^Z{t,x)X{R<l}(t)\\LP.<^ < t““, 


which finishes the proof. 


□ 
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6.1. Proof of Theorem 12.351 Now we are ready to prove Theorem 12.351 

Proof of Theorem \2.S5\ We split the proof into two parts. We first study the esti¬ 
mates for Z. The estimate for Y is substantially more involved and we do it after 
studying Z. 

The estimates for Z: The strategy of the proof here is the same as in [sni p. 14, 
15]. Suppose first that uq G is such that |a;| |“o(2;)| dx < oo. By Lemma 

13.41 there exists f) G such that 

Mo = Minitdo + davf) 

and \\(j)\\L^ < C'd|||a;|Mo||Li- Consequently, 

u^mt{t,x) = Minit {Z{t, ■)-kdo) {x) + {Z{t, ■) * div^(-)) {x) 

= MinitZ{t, x) + (yzit, ■) -k 4>){x), 

which yields 

Uinitft, x) - MinitZit, x) = iyzit, •) * 4>){x). (6.26) 

By Young’s inequality it follows that for any 1 < p < ki(/3, d) 

< ||VZ(f,.)||Lp||<^||l < ||VZ(t,.)||Lp|||x|uo||Li 

where we used Lemma [6.221 Hence 

t^ •) - MinitZ{t, •)||lp < ft, 

which is the first part of assertion (ii). The second part follows from (16.2611 by 
applying Young’s inequality for weak L^’-spaces [231 Theorem 1.2.13]. 

To prove (i) we choose a sequence {r]j) C such that Pj dx = Minit for 

all j and rjj uq in L^(]R^). For each j by Part (a) and by Lemma I^T] we obtain 

*) MinitZ(t^ ‘)I|ap 

< \\Z{t, •) * (uo - Vj)\\LP + \\z{t, ■)-kr]j - MinitZft, Ollip 

< \\z{t, oiUpIIuo-p.IIli + c{j)t-f-t (1-1) 

< Cit~~ (^“p)||uo - Pjllii + C{j)t~~P~~F (i-p), 

and therefore 

i~(l-p)||Ui„it(f, •) - MinitZ{t,-)\\LP < Clllwo - Pjllii +C[j)t~'^, 
which implies 

limsup •) - MinitZ{t, •)||lp < Ci||mo - PjIIli- 

t—^oo 

Assertion (i) follows by sending j ^ oo. This finishes the decay estimates for Z. 
We continue with Y. 


The estimate for Y: Next we turn to study Uforc- We split Mforc into two parts 
as follows 


Mforc 



f{T,y) dy dr + 


■OO 


f{T,y)dydT 
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and note that 




Y{tr) [ [ f{'r,y)dyd- 

Jt JR'i 

poo p 

/ / |/(r,y)|dyd 

Jt Jr<^ 

/ \f{T>y)\dydT ^ 0 


as t —>■ cxD. Here we used Lemma IB. 121 to obtain 

\\Y{t,-)\\L. 


t^oo. 


Therefore it suffices to prove that 


[ iY{,t-T,-)-k f{T,-))dT [ f f{T,y)dydT 

Jo Jo Jr'^ 


Lp 


as t 


To prove the assertion, we fix 0 < <5 < and decompose the set of integration 
(0, t) X into two parts 

n,it) = iO,St) X {y : |2/| < (5^/^}, 

^2(0 = ( 0 , t) X \ 


Let us start with the set We estimate by using the integral form of the 

Minkowsky inequality in the case 1 < p < 00 to obtain 



-T,--y)-Yit, •)]/(t, y) dy dr 


-T,- -y)-Y{t,-) 


\f{T,y)\dydT. 

Lp 


(6.27) 


If p = 00 , the same estimate holds trivially. Note that in Hi(t) we have t > t — t > 
t{l — S) > so t — T and t are comparable and there is no singularity in t. Our 
aim is to prove that the L^’-norm on the left-hand side of (16.2711 tends to 0 as d ^ 0 
uniformly in t. To achieve this, we distinguish two different cases w.r.t. x € 
when looking at the L^’-norm on the right-hand side of (16.271) : 


(i) \x-y\< 

(ii) \x — y\ > 2(i5t)“/^. 


Observe that this splitting seems to be needed. If we simply estimate the L^’-norm 
by the triangle inequality 


Y{t-T,--y)-Y{t,-) 
we would get a bound 


LP 


LP 


< Y{t-T,--y) 


Y{t,-) =:/i+/2 , (6.28) 

Lp 


which is of a right form but the problem is that this quantity does not converge to 
zero as i5 —?► 0 which is what we are after. Therefore we need to do the estimates 
more carefully. 


The motivation for the splitting is that in the case (i) both \x — y\ and |a;| are 
bounded from above by a multiple of ((5t)“/^. In this case we will simply use the 
triangle inequality (I6.28p . The second case (ii) is more complicated, but here we 
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will proceed as follows. Since there are differences both in the space and the time 
variable, we treat the differences separately by using the triangle inequality: 


r(t-T, • -y) -y(t, •) 


< 


Lp 


Y{t-T,- -y) -Y(t-T,-) 
Y{t-T,-)-Y(t,-) 


Lp 


+ 

=: h + h- 


(6.29) 


LP 


In both of these we shall use the Mean Value Theorem. Note that in the case (ii) 
we are away from the singularities of x and f, since t — T> and |x| = \x — y + y\ > 
\x-y\- \y\ > 

Since the asymptotic behavior of Y is different for each d, we consider here only 
the case d > 2(3 and 0 < a < 1. The other cases can be treated similarly, since 
the proof is based only on the pointwise estimates for F, VF, and dtY given in 
Lemmas KE\ and K7\ 

We start with the case (i). Note that for (r, y) G fli(t) we have 

\x — y\ 2(5“/^ 

[t _ tY/P - (1 - (5)“//3 ’ 

so we may use the asymptotic behavior of Lemma 14.61 for small values of the simi¬ 
larity variable R to obtain 

\Y{t,x)\<t-‘^-^\x\-<^+^^. (6.30) 

As mentioned before, we use (16.2811 and (16.301) to obtain 

^ J\x-v\< 2 {&t)P‘/P ' 

Introducing the spherical coordinates gives the desired estimate 

/i < 

Notice that the assumption p € [1, K 2 ) guarantees the integrability and the positiv¬ 
ity of the power of 6 , which is needed in the end. The same proof applies also for 

h- 

Now we shall provide the estimate in the second case (ii). Since we are going to use 
the Mean Value Theorem, we need to calculate the derivatives of the fundamental 
solution F. We recall the following estimates from Lemma [4.71 for d > 2/3: 

|VF(t,x)| 

and 

\VY{t,x)\<t-°‘-^\x\ 

By using the Mean Value Theorem for I 3 we obtain 

h = |y|||VF(t-T,i(-))||Lp 

for some x on the line between x — y and x, where x denotes the integration variable. 
Since 

\x\ = \x - y + X - {x - y)\> \x - y\ - \x - {x - y)\ 

>\x-y\- \y\ > 


> 1 , 

(6.31) 

< 1. 

(6.32) 


(t-T-y/P - 2t°‘/P 




we have 


(6.33) 
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Notice, that since S can be small, we have to use the asymptotics near zero and near 
infinity. Therefore we divide the integral I 3 into two parts /31 and I 32 depending 
on whether — t)““ is less than 1 or greater than 1. 

In I 32 we use 

3 

|i| < \x-y\ + \y\ < ^\x-y\, 
so the set of integration is contained in the set 

{xeM'" : \x-y\>^{t-Tr/^}, 

which implies the estimate 


I 32 < ( [ {t- dx] 


1/p 


'\x-y\>l{t-T) 
< P+2a-l 


\x-y\>l(t-T)'=‘/P 


\x - y\^-^-^-P'>P . 


Introducing spherical coordinates gives the estimate 

which is of the form we need. 


For J 31 we note that by (16.3311 the set of integration is contained in the set 
(a; G < l}, 

I - 2 {t- - J 

so by using (16.321) we obtain 

I / 




/(,< \x-v\ 

- 2{t—r)^/^ - 


Once again we use the fact that |i| and \x — y\ are comparable. We may proceed 
as before except we have to separate two cases: (a) {—d — I + 2l3)p = —d or (b) 
{—d— 1 + 2/3)p ^ —d. An easy calculation shows that the first case is possible only 
in the case /I > ^. The case (a) leads to a logarithmic function. Indeed, we may 
estimate 


/31 < 

< (5“/^|log(5r''V/^"““i 


- 2 (t-T)“/^ - 


\x - y\i-d-l + mp dx^ 


1/p 


A simple arithmetic calculation shows that the power of t is actually 

a ad 1 

which is exactly of the right form and the factor depending on S tends to zero as 
(5 —>• 0 uniformly in t. 


The assumption p G [1, ^ 2 ) leads to a usual power function similarly as before. We 
omit the details and write the final estimate 


-^31 ^ 


5^ -5^ 




(6.34) 


Again the assumption p G [1, ^ 2 ) guarantees that the second power of 5 is positive, 
so we have obtained the desired estimate also in this case. 
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For /4 we use again the Mean Value Theorem to obtain 

h = T\\dtY{i,-)\\L. 

for some t € {t — T,t). Note that in r2i(t) t and t are comparable: (1 — 5)t < t < t. 
Now \x\ = \x — y + y\ > \x — y \ — \y\ > so 


= ^«//3 


fa //3 — la/p — (6.35) 

and again we have two cases, since 5 can be small. We denote the integrals by /41 
and I 42 depending on whether 5 < 1 or 5 > 1 . 

Again, we recall from Lemma 14.71 the estimates 


\dtY(t,x)\<t-^-^\x\-'^+^^ ^<1 


1 ^ 


and 




\x\^ 


> 1 . 


(6.36) 

(6.37) 


The estimates (16.3511 and (|6.36l) now give for /41 that 


/41 ^ 8t 

By changing the variables x -fA 

hi < 

< 


toe/13 > |x| >((5t)“/^ 

=: z, we obtain 


ih°‘-^)p\x\hd+mp 


( [ \z\hd+mPdz) 


i/p 


ijp 


/ 5 “/' 3 <|z|<l 

r~~ 




Since the powers of 5 are even better than in (I6.34L we have derived the desired 
estimate for / 41 . 

For J 42 we observe that 

1*1 

1 < ; <_!_!_ 

- - ((l-, 5 )t)“/^ 

which implies 

\x\>{{i-5)tr/h 

We use (16.371) to obtain 


h2 'h St 


1^1 


^'(2a-2)p|^|(-d-/3)p^^\ 

/ 


1/p 


Making the obvious change of variables x •(-)• =: z we end up with the estimate 

/42 


similarly as before. 

Collecting all above we see that 


,l+^{l-i)-a 


II ' 


{Y{t-T,- -y)-Y (t, •))/(*, y) dy dr 


< 


LP 


S’^Wfh 


for some positive number rj. The upper bound tends to zero as 5 ^ 0 uniformly in 
t. 
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We now fix (5o < 5 such that the previous term is small and continue to estimate 
the norm 




^2{t) 


{Y{t-T,- -y)-Yit, ■))f{T, y) dy dr 


LP 


Using the integral form of the Minkowsky inequality we have 




(Y{t - T, ■ - y) - Y (t, •))/¥, y) dy dr 

p^~°‘ [f Y{t-T,--y) \f{T,y)\dydT 

LP 


='h + h 




n 2 (i) 

n 2 (t) 

( 

02 (i) 


LP 


yit,-) |/(r,y)|dydT 
LP 


By Lemma [6.121 we have that||y(t, •)||lp ^ p^ ^ and, therefore, we may 

directly estimate Iq by 


as t ^ oo. 


h < 



\f{T,y)\dydT 0, 


For J 5 we have two possibilities: either r < Sot or r > Sot. According to this we 
split the domain D,^\t) into two parts: 

uf (<) = {0,Sot) X { 2 / e R'" : \y\ > U (doi, <) x R'', 


where ( 0 ) indicates the fact that we have fixed d = i5o. 
Hence, can be written as 


A =tl+¥(l-F)-- 




f'Sot 




||r(t - r, • - 2 /)||lp|/(t, y)\dydT 


[ f -y)\\Lp\fiT,y)\dydT. 

J 6 ot 


We use the same bound ||y(t, •)||lp < ^ as above for both integrals. 

Then the first integral is dominated by 

f f (t - T)°‘~~^^^~py^\f{T,y)\ dr dy 

Jo J\v\>{Sot)P‘/P 

<il-So)°‘~'^^^~py^ f [ \f{T,y)\dTdy, 

Jo J\y\>(5ot)'-/^ 

which clearly tends to zero as < —> oo. The upper bound for the second integral is 

J Sot 


This integral causes problems, since now there is a singularity in t. But the as¬ 
sumption p S [1,R2) guarantees that the singularity is weak. We use the decay 
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condition (12.341) imposed for the source term. By using this, we have 

J ^ 0 * 

p^Qt 

Jo 

which tends to zero as t —?> oo, since 7 > 1. This, finally, finishes the proof of the 
case where d > 2/3 and 0 < a < 1. The other cases are proved similarly. We omit 
the details. □ 


7. Optimal T^-decay for mild solutions 


In this section we will give the proof of Theorem l2.36l Here we only consider equa¬ 
tion ( 1 ^ . but our reasoning can be extended to cover a wider range of equations. 
The main tool we use is Plancherel’s theorem, but in general it can be replaced by 
more general multiplier theorems which allow one to study more general equations, 
too. For details of such an approach we refer to our earlier paper |28j . Here we 
restrict our study to equation (12.51) for simplified exposition. 

We begin this section by showing that our decay rate is optimal. Indeed, we have 
the following result. 

Proposition 7.1. Let a G (0,1), d > I, and d ^ 2/3. Suppose u is the mild solution 
of the Cauchy problem ( 1 ^ with / = 0. Assume further that uq G nT^(R'^) 

with uq da; 7 ^ 0. Then 

t>l. 

The eonstant in the estimate depends on uq dx. 


Proof. Let po > 0, t > 0 and p = p{t) G (0, po]- By Plancherel’s Theorem, mono¬ 
tonicity of and the estimate Eap{—x) > ci/(I -I- a;) for all a: > 0 (with some 
Cl >0), we have 


\\u{tr)\\l.= 


> 


> 


\\uit,.)\\l.= / \Z{t,0f\M0\^d^ 


(27r)“ 


'Bp(O) 


Eo.A-\^fn^\MOfdf 


(27r)'^(I-h p/3t“)2 

C2 




(1 -bp^'t“)2 


'J Bn 


(7.2) 


By the Plancherel Theorem and the Riemann-Lebesgue Lemma we have uq G 
Co(R‘^) n L^(]R‘^). By the Lebesgue differentiation theorem, we may choose po 
small enough in order to obtain 

P~‘^ [ |mo(C)P^C > for all p e (0,po]. 

Using this in (17.21) gives the lower bound 


\Htr)\\h> 


C2|u(0)pp^ 

2(1-h p/5t“)2 ■ 


(7.3) 
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Next we choose p = po, which yields 

for t > 1. On the other hand, the choice p = p{t) = gives < p^ 

and thus by (IZ31) we get the estimate 

t>i. 

These estimates combined together give the claimed lower bound. □ 


Observe that the constant in the above proposition is of the form C = C (po) I fgd udx\, 
where also po depends on | udx\. Nevertheless, we obtain that the decay rate 
in Theorem 12.361 is optimal. We will now give a proof of this decay result. 


Proof of Theorem \2.36l To prove the upper bound, we proceed as in Theorem 
4.2]. Suppose that d < 2/3. By Plancherel’s Theorem, the Riemann-Lebesgue 
Lemma and the estimate (|2.22|1 . we have 


Ht,-)\\h = Mtr)\\h= / \Z{t,0\^\uo{0\^d^<\\uoU 


f \zit,0? 

Jb.'^ 


< ll^oll 


dr] 


(i + iei/5t“)2 


/R'i 


(1 + |#) 2 ’ 


d^ 


(7.4) 


where in the last step we have made the change of variables ^ -fA =: p. Now 
the condition d < 2/3 guarantees that the last integral is converging. Hence we have 
derived the upper bound in the case d <2/3. 

We are left with the case d > 2/3. Here we use the Hardy-Littlewood-Sobolev 
Theorem on fractional integration. Indeed, we choose g = 2 in Theorem 13.51 to 
obtain 

||(-A)"^uo||l2 < IImoII^^ < oo, (7.5) 

since uq S L^(R'^) fl L^(R'^) implies that uq G L'3+?(R'^) by interpolation. Using 
this and the estimate 12.211 we have 


Mtr)\\h = f \e^\z{t,om\-^Mo\^d^ 

Jk.’‘ 

Jr’‘ 


which completes the proof by (EH). 

For the borderline case d = 2/3 we estimate directly by Young’s inequality (13.31) to 
obtain 


||u(t,.)||L 2 .=. = \\uo{-)*Z{t,-)\\L 2 .^<C\\Z{t,-)\\L^.^\\uoh^<C\\uoh^t-^, 


where we used Lemma IQ to estimate the weak L^-norm of Z. This finishes the 
proof. □ 
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8. Energy method and L^-decay for weak solutions 

In this section we consider the L^-decay of weak solutions which are defined in 
Definition 12.91 We will restrict our study to the homogeneous case / = 0. We will 
proceed in a rather formal manner where we prove the estimates starting directly 
from the equation by multiplying it with the appropriate test functions. For the 
details required for the rigorous treatment starting from the Definition l2.91 we refer 

to [55]. 

In the proof of Theorem 12.371 we will need the following Lemma from [43] . 

Lemma 8.1. Let T > 0 and LI C be an open set. Let k G kFioc ([0, oo)) be 

nonnegative and nonincreasing. Then for any v G L^((0,T) x 12) and any iiq G 
Lf{Ll) there holds 

[ vdt{k * [v - Vo]) dx > |k(t, •)llL 2 (a)at(fc* [||i;||l 2 (q) - ||i;Q||i 2 (n)])(t), (8.2) 

Jn 

for almost every t G (0,r). 

Proof. The result is originally from [43] . For the proof in our context we refer to 
Lemma 6.2 in [55]. □ 


Observe that in our case the kernel k corresponds to gi-a- The function gi-a is, 
however, not in W^’^. For this reason, a rigorous treatment of the problem requires 
an appropriate regularization of the fractional derivation operator in time. One way 
to do this is via its Yosida approximations, which leads to an integro-differential 
operator of the same form with a kernel gi-a,n that is also nonnegative and nonin¬ 
creasing, and which belongs to oo)). The details of such calculations can 

be found in [55], see also [25]. Note that the regularized weak formulation used 
in [28] and [46] does not involve an integral in time on [0,r], but it requires the 
validity of a certain relation pointwise a.e. in (0,r). Here we proceed on a formal 
level by using the singular kernel gi-a and a formulation of the problem where we 
only integrate in space (not in time) against a test function. 


Lemma 8.3. Let uq G L^(R‘^) n L^(R'^). Suppose u is a weak solution of equa¬ 
tion (12.81) with initial condition M|t=o = uq, and assume that \2.38\) is satisfied. 
Then 


for a.a. t > 0. 


\\u{t, •)IIli(R‘^) < llnollii 


Proof. Letting i? > 0 we choose a nonnegative cut-off function if G CqIBh+i) such 
that = 1 in Bn and i/j < I as well as IVi/jI < 2 in Bn+i- Here Bp denotes the 
ball of radius p > 0 and center 0. For e > 0, define 

He{y) = (.y^+ 2/e K. 


Clearly G C^{ 


and Hf G 

H'M = 


Indeed, 

y 


{y^ 


yGR. 


Observe that is convex. Testing the PDF with Hf{u)if gives 


/R<1 


Hf{u)lfd^{u — Uq) dx + Feit) = 0, 


t > 0, 
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where 

Fe{t)= [ [ K{x,y)[u{t,x) - u{t,y)][Hl{u{t,x))ip{x) - Hl,{u{t,y))ilj{y)] dxdy. 

JR ‘1 Jmd. 

Since Hg is convex, we may (formally) use the inequality from Corollary 6.1 in |28j . 
to the result that pointwise a.e. we have 

H'^iu)d?{u - uo) > d?{H,{u) - H,iuo)). 

Applying this to the previous relation and convolving the resulting inequality with 
ga we obtain 



{Hs{u) - Hs{uo))'tpdx + ga * <0, 


t > 0. 


Next, we send £ —?► 0 and observe that Hg{y) \y\ as well as H'^{y) —>• signy for 
y GR. Thus we get 


x)\ - \uQ{x)\)ii{x) dx + {ga * F){t) <0, t > 0, 


(8.4) 


with 


m 


f / 

Jr** JR<i 


K{x,y)[u{t,x) 


-u{t,y)] ■ [signM(t,x)^(x) - sign u{t,y)'ijj{y)]dx dy 


> / K{x,y)[u{t,x)-u{t,y)]-[ip{x)-ip{y)]signu{t,y)dxdyFi{t). 

7r<j jRd - 


Using the properties of ^ and K, we may estimate as follows. 

f \u{t, x + h) — u{t, a;)| ■ \ip{x + h) — ip{x)\ 


l^iWI <A 


= A 


< 


/R ‘1 jR'i 

I f 

JR'i J\h\>R/2 

2A f f 


|/j|d+/3 


.. .dhdx + A 


dh dx 


.. .dhdx 


jRd\BB/2 J\h\<R/2 

\u{t, X + h) — u{t, x)| 


(i?/2)2 J\h\>R/2 


\h\ 




dh dx 


+ A 


\u{t,x + h) — u{t,x)\■ 2^ 2 |Vi/’| 


/RASr/2 Ah\<R/2 


\h\ 


d+i 


■ dh dx 


< 


2 ^+i A 


[u(t, •)! /3 + 2A / 

^ ^ ^V 2 a(R.) 


\u{t, x + h) — u{t, a;)| 


R‘^\Br /2 -tR-^ 


\h\ 


d+i 


dh dx, 


where the last two terms tend to 0 for a.a. f > 0 as i? —>■ oo, by assumption 
(12.381) . Thus the assertion follows from (18.41) and the previous estimates by sending 
R —>■ oo. □ 


Finally, we have the following Lemma, which shows that the L^-norm of a weak 
solution is a subsolution to a purely time-fractional equation. 

Lemma 8.5. Let K, uq, and u he as in the previous lemma. Then there exists a 
constant p = /i(d,/3, A, ||uo||li) > 0 such that (formally) 

9?'[\\u{t,-)\\L2 -\\uo\\l 2 ] + p\\u{t,-)\\j(t^ <0, t>0. (8.6) 

Proof. Choose the test function ip = u in Definition l2.9l of weak solutions. We apply 
Lemma 18.11 for the fractional time derivative to obtain 

ll^*(^c)IU25(“[||u(t, •)|U2 - ||mo||l=] + [ [ K{x,y)[u{t,x)-u{t,y)]'^dxdy <0. 
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For the elliptic term, we then use the fractional Nash inequality (cf. [HI p. 13]) 
together with the assumption (1^ on the kernel K and with Lemma 18.31 to obtain 


< C'||u(t,-)|| 


20 

-T 


20 

< C\\uo\\J, 


I I I \d+B 

's.'i JRd' \x — y\‘‘-+P 

/ K{x,y)[u{t,x) - u{t,y)]'^ dxdy. 


This concludes the proof. 


□ 


8.1. Proof of Theorem 12.371 We are finally ready to prove the decay result for 
the weak solutions. The proof is based on using the comparison principle for the 
purely time-fractional equation (18.611 . 

Proof of Theorem \2.37\ Let T > 0 be an arbitrary real number. By the comparison 
principle for time-fractional differential equations (see [431 Lemma 2.6 and Remark 
2.1]), the inequality (18.6|) implies that \\u{t, •)||l 2 < w{t) for almost every t G (0, T), 
where w solves the equation corresponding to (18.6|) . that is 

dfiw - WQ){t) + y,w{t0 = 0, t>0, w(0) = Wo := ||mo||l2, 

where we put 7 = 1-1-^. It is known that for wq > 0 there exist constants 
Cl, C2 > 0 such that 

Cl , , C2 

- W < < -IT: t>0, 

l-|-t^ l-l-t^ 

see [321 Theorem 7.1]. Since T > 0 was arbitrary, we conclude that 

C2 C2 

||w(t, < w{t) < -^ almost every t > 0. 

I + t~i 1 + t d+ys 

This finishes the proof of Theorem 12.371 □ 
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